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We present the study of parametric resonance in a one-dimensional cavity based on the analysis of 
classical optical paths. The recursive formulas for field energy are given. We separate the mechanism 
of particle production and the resonance amplification of radiation. The production of photons is a 
purely quantum effect described in terms of quantum anomalies in recursive formulas. The resonance 
enhancement is a classical phenomenon of focusing and amplifying beams of photons due to Doppler 
effect. 
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^ ' I. INTRODUCTION 

o ' 

QThe phenomenon of parametric resonance is related to the instability of open systems under the action of an external 
periodic force. In quantum field theory, we expect the resonant amplification of quantum fluctuations 0, 0, Q • The 
resonance enhancement of vacuum fluctuations is usually referred as the dynamical Casimir effect. The standard 
C*~) ■ model to investigate this effect is the system of electromagnetic fields confined inside a vibrating one-dimensional 
cavity ilia. 

Meplan and Gignoux [8j found the correspondence between the wave propagation in a vibrating cavity and the 
motion of massless particles in a two-dimensional space-time billiard. They used the generalized Korringa-Kohn- 
Rostoker method. The method presented in this paper offers both a better insight in resonance mechanism and a 
• more useful tool for detailed calculations. Moreover, the brief discussion of the quantum version in § does not trace 
£NJ \ properly the impact of conformal anomaly [jj . 

In this Letter, we present a new approach to study resonance solutions. We show that the resonance amplification 
together with the formation of narrow packets in energy densities [Tel ] is a purely classical phenomenon. It can be 
explained as a consequence of a cumulative Doppler effect . In this way, we are allowed to interpret a wave packet 
as a beam of massless and non-interacting particles. Hence, we introduce the billiard function which contains full 
Qh, information about possible reflections from a cavity wall. This function is helpful to establish recursive formulas for 
■ physical quantities in a vibrating cavity. Production of particles in quantum case is represented by conformal anomaly 
contribution in such formulas for energy densities. Our technique is illustrated with solutions of several examples of 
vibrating cavity systems discussed in literature. 

a*' 



II. CLASSICAL BILLIARDS 



Before studying the field-theoretic models, we consider briefly the classical billiards. Let us start with the simplest 



example, namely a head-on collision of two non-relativistic objects. A particle moves with velocities v and v' before 
and after the collision, while respective velocities of a target are denoted by u and v! . This elementary physical 
problem is fully described by the reflection law in one dimension (on a line), 

v + v' = u + u' . (1) 

If we assume a target to be very heavy and skip its recoil, then the velocity of a particle after the collision is just 
v' = —v + 2u. We sec that if a particle were bounced with a regular frequency, then its energy would grow quadratically 
with time. Let us look now at a head-on collision of relativistic objects. This two-body problem is described by the 
relativistic reflection law which says that the sums of rapidities are equal, 

artanhu + artanhv' = artanhu + artanh it' . (2) 

In this case, if a particle bounces rhythmically on a heavy target, then its rapidity grows linearly with time. It follows 
that the energy grows exponentially with time. In particular, we can consider a left-moving massless particle (photon) 
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Figure 1: The definition of the billiard function /(«). It returns full information about possible reflections of a light-like 
particle from the mirror moving with a prescribed trajectory L(t). The billiard function identifies the trajectory uniquely and 
its derivative is the retarded Doppler factor. 



hitting a target with a velocity of light. The exact formula for the change of the particle energy after a single collision 
yields: 



1-u / 2E l-u\ 

E' = E 1 + — \ . 3 

1 + u \ M V 1 + u J y ' 

Two consecutive factors on the right-hand side are obviously the Doppler factor and the Compton factor. For a large 
target mass M, the Compton factor is neglected. Therefore, during each head-on collision the energy of a massless 
particle increases (or decreases if the target moves in the same direction) by the Doppler factor. 

Now, we insert a heavy target (mirror) which moves with some prescribed trajectory L(t). Let us then introduce 
the following billiard function\\'^ : 

f(t + L(t))=t-L(t) . (4) 

We define the retarded time t* using the retardation relation t = t* + L(t*) and recognize that the derivative of the 
billiard function is just the retarded Doppler factor: 

/(*) = ■ (5) 

W 1 + L(i*) V ; 

Since the derivative is additive, the billiard function is increasing and its inverse is well defined. The Doppler factor 
is greater than one (amplifying) if the mirror is moving left, i.e. towards an arriving test massless particle. For any 
physical trajectory of the mirror, its corresponding billiard function is denned unambiguously. On the other hand, any 
proper billiard function may be used to reconstruct its genuine mirror's trajectory. More information can be found 
in |l2( . Although rather cumbersome in general, the billiard function is the most useful characteristic to analyze the 
parametric resonance within a particular cavity model. We can make use of the correspondence between waves and 
massless particles and relate the field theory problem with the analysis of trajectories in a space-time billiard @- 



III. OPTICAL CAVITIES 



Under the guidance of the above simple examples from classical billiards, we consider the case of field theory. As 
usual, we will discuss the ideal resonator model [4| formed by two perfect mirrors. The electromagnetic potential 
A(x, t) obeys the one-dimensional wave equation and any classical solution is composed of left- and right-moving wave 
packets, 

A(t,x) = A L (t + x) + A R (t-x) . (6) 

For the sake of simplicity, we assume the left mirror to be fixed at x = 0, whereas the right mirror is vibrating according 
to a prescribed trajectory L(t). Then, the field A(x,t) is subject to Dirichlct's boundary conditions, namely: 

Al(t) = —Ar(t) = ip(t) , <p{t)=<p(J{t)). (7) 

We follow the solution of the Cauchy problem with the help of the billiard function. From initial data for A and dt A 
at t — 0, we read the function tp for arguments in the interval [— L(0),L(0)] (up to some irrelevant constant). The 
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solution is obtained if we extend ip into the whole real domain. We can do it using Eq. J7J provided that the billiard 
function is known. Usually, if we predetermine some trajectory of the mirror, then the corresponding billiard function 
is to be found only in the way of numerical iterations based on Eq. (@J. On the other hand, to specify the mirror 
motion we can predetermine the billiard function instead. 
The energy density of the classical wave packet is given by: 

T 00 (t, x) = l/2(8 t A) 2 + l/2(d x A) 2 = g(t + x) + g(t - x) , (8) 

where q(t) = ip 2 (T). The formula for the total energy can be presented in the following form: 

fL(t) rt+L{t) 

E(t) = / dx T 00 (t,x) = / dr q(t) . (9) 

JO Jt-L(t) 

The billiard function determines all possible trajectories of a massless particle or all classical optical paths inside the 
cavity. Each trajectory can be represented by the sequence {T„(r)} = {(/ _1 )°"( T )} n=0 , where (/ _1 )°™(t) denotes 
here n-fold composition / _1 o / _1 o . . . o / . A light-like particle starts at time r moving right from x = and next 
elements of the sequence read next times of collisions with the static mirror. We define also retarded times T*(t) by 
demanding T*(r) + L(T*(r)) = T„(r). It is easy to find the relation T*{t) = (T„(r) + T„_i(t))/2. 

The appearance of the parametric resonance is related to the existence of periodic particle trajectories 0- Such 
trajectories obey the following periodicity condition (for any non- negative integer n): 

T n (T ) = T + nT , (10) 

where tq is a starting point in time and T is a period. The retarded times are respectively T*(tq) = T„(to) — T/2. 
Looking at Fig.l, we notice immediately that any periodic particle trajectory appears on condition that the trajectory 
of the mirror has return points [l3T ] : 

L (tq + nT) — T/2 . (11) 

It is usually assumed that the position of the mirror T/2 refers to the length of the static (unperturbed) cavity, so 
that T/2 = L = L(0) in this paper. 

Let us define the characteristic for particle trajectories which will play the most important role in our analysis. It 
is the cumulative Doppler factor: 

D n {r) = J— = ft /(T„(t)) = f[ 1 ~ L ™ T)) . (12) 
Tn{r) t=\ fJlHKW) 

We will call a trajectory to be positive (stable, attractive) if its cumulative Doppler factor tends to infinity for large 
n, and to be negative (unstable, chaotic) if its cumulative Doppler factor goes to zero with increasing n. 

The existence of periodic particle trajectories inside the cavity is equivalent to the existence of suitable return points 
in the trajectory of the cavity wall. In particular, a periodic particle trajectory is positive (negative) if any return 
point of the mirror trajectory gives a Doppler factor greater (less) than one. It happens if the right mirror is always 
moving toward (outward) the left mirror at all return points. 

Denote by r + and r_ starting points for positive and negative periodic trajectories respectively. For small pertur- 
bations e, it is straightforward to derive the following formulas: 

T„(r+ + e) = t+ + nT + eD- 1 ^) , (13) 



T n {T^+eD n {T-))=T-+nT + e . (14) 

Now, we can show how the existence of periodic particle trajectories inside the cavity triggers the resonance 
instability of the system. For the sake of simplicity, we assume that there are only two periodic particle trajectories, 
a positive one and a negative one, with starting points r+ and t_ lying in the initial interval [— L(0), L(0)]. These 
simplifications are only for clarity, and the method can be adapted for more sophisticated examples. As it was said, 
to solve the Cauchy problem we need to extend the function ip given in Eq. (JJJ from the initial interval to the whole 
real domain. If we are interested only in the evolution forward in time, then it is enough to obtain data only for 
arguments to the right of the initial interval r > L(0). To build such extension of tp, we pick points from the initial 
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interval, run along particle trajectories and read extended values of tp from Eq. J7J. To analyze our solution obtained 
in this way, let us look first at the profile function of the energy density. If we take the following iteration formula: 

e(T„(r)) = Q(T)D 2 n (r) , (15) 

together with Eqs. (|13I14|) . then it is straightforward to derive the following asymptotic formulas for large arguments 
(n > 1): 

e (nT + t++ eD-\t+)) S q{t+ + e)D 2 n (r+ + e) , (16) 



g{nT + r_ + e) = q{t- + eD n { T -))D 2 n (r- + eD n (r^)) . (17) 

The above formulas explain the formation of travelling narrow packets in the energy density T 00 (t, x) 0, 0, ITc| . 
The profile function q(t) is concentrated around spots of the positive periodic trajectory. From Eqs. H16I17JI one 
can estimate easily the width and the height of peaks. The height grows like D 2 1 (t + ) : and the width diminishes like 
D~ 1 (t+). The values of the energy density far off the peaks decreases according to D^(r_). 
The total energy can be analyzed using the following formula: 

pT„(to) pro 

E(T*(t )) = / dr q{t) = / dr Q{r)D n {r) . (18) 

JT„„ 1 (To) Jf(T ) 

We can calculate the evolution of the total accumulated energy using only initial data and the cumulative Doppler 
factor. It is apparent that an initial shape of the classical field inside the cavity is of minor importance. Let us test 
our method with an example of sinusoidal cavity wall motion: 

L(t) =L + Aisin(wi) . (19) 

Obviously, we are to assume that AL < L and ojAL < 1. The parametric resonance frequencies are ojn = Ntv/L, 
where N is the order of the resonance. It is easy to find periodic particle trajectories. First, we read all return points 
of the mirror from Eq. ljTT|l : sin(wjvT*) = 0. They correspond to the following starting points of positive and negative 
particle trajectories: 

_ (-N + 2m+l)L _ (-N + 2m)L 

where m — 0, 1, N — 1. The corresponding values of the cumulative Doppler factors can be calculated: 

D n (T± m ) = ( ] ±UJN ^ } = exp (±2nuj N AL) for AL <C L . (21) 
V 1 T lo n ALJ 

There are N travelling peaks in the energy density. All the peaks have the same height and the same width. The 
corresponding time evolution of peaks agrees with the results of quantized version [6j. The formation of peaks is 
basically a classical phenomenon. 

The resonance instability in the field theory appears not only for finely tuned frequencies. We obtain usually some 
band structure. Our method enable us to get insight into off resonant behavior of the cavity system as well. Therefore, 
we assume some perturbation of the resonant frequency: ui — ujn + Alj. The equation for return points is now: 

sin (l,t*) = — . (22 

uiAL 

The solutions for return points exist provided that: 

— <—. 23 

to L 

This condition has been obtained numerically for quantized version in |14| . It defines the band structure of the 
parametric resonance |15|. For small perturbations, the structure defined by Eq. ll'.'il is typical and resembles those 
obtained in classical mechanics. The cumulative Doppler factor at starting points yields: 



HvPIFW (24) 
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For off resonant vibrations of the cavity, the formation of narrow peaks and the energy growth take a longer time. 
Let us conclude our classical analysis. Within the framework of our method the calculations are straightforward 
for any type of cavity motion. The parametric resonance is related to the existence of periodic trajectories in a 
two-dimensional cavity billiard. As it was already noticed in |13|. there is no need for periodicity of cavity motion. 
We need only periodic returns of the cavity to the unperturbed position. For typical cavity motions, the growth 
of the total energy is exponential with time. The cumulative Doppler factor describes all details of the formation 
of narrow wave packets and the energy growth. Essentially, many details of the resonant enhancement of classical 
electromagnetic radiation inside a vibrating cavity match the results of quantum theory. 

Finally, we discuss vibrating optical cavities in the quantum field theory. Apart from the interference effects which 
we have already demonstrated in the classical theory of the cavity, we are to account for instability of the quantum 
vacuum and possibility for particle production. The vacuum expectation value of the energy density is given by 3 : 

(T 00 (t,x)) = g(t + x) + g(t-x) , (25) 

where 

e(r) = ^ 2(r) ~i 5[i?](T) • (26) 

The second term is responsible for particle production. It is defined by the Schwartz derivative: 

sm = 5-5(5)' • (27 > 

The crucial information is contained in the phase function, which obeys the following Moore's equation (being a 
counterpart of kinematical laws Eqs.(l,2[7J): 

R(r) - R(f(r)) = 2 . (28) 

We see that our billiard function is useful for the quantum case as well. It is straightforward to derive the following 
recursive relation: 

g(T n (r)) = q(t)DI(t) + A n (r)D 2 n (r) , (29) 

where the cumulative Doppler factor is given again by Eq. (|12fl and the cumulative conformal anomaly contribution 
can be calculated from the following formulae: 

1 1 " 

An{T) = 2^ S[Tn]{T} = -^E^ 2 ( r ) 5 [/](r fe (r)) . (30) 

fc=i 

The total energy can be calculated now from: 

E(T*(r )) = / " ( ° } dr g(r) = f° dr [g(r) + A n (r)} D n (r) . (31) 

•/T n _i(r ) Jf(r ) 

In comparison with the classical formula Eq. I|15|l . aside from the Doppler factor there appears a new additive term 
on the right hand side. It represents the conformal anomaly of the theory. Most of classical results, like the band 
structure around resonance frequencies, the formation and the shape of travelling packets in the energy density, the 
exponential growth of the total accumulated energy, can be reproduced in the quantum case as well. These features 
are to be calculated from the cumulative Doppler factor. The novelty of the quantum description is the conformal 
anomaly. We show that for the lowest resonance channel, the anomalous mechanism of energy growth clashes with 
the resonance enhancement of the initial vacuum fluctuations. 

Let us refer again to the sinusoidal motion of the cavity l|19|) . We consider motions with resonance frequencies wjv- 
At the beginning, the cavity is static and empty. It corresponds to the following condition: 

2 

q(t) = ^ = -- ^- for re [-L, L] . (32) 

^ V ' 48L 2 48tt L ' J v ) 

The above value is just the static Casimir energy density which is present even the cavity is at rest. We can compute 
the cumulated anomaly contribution at starting points of positive trajectories directly from Eq. I|3U|) : 

Mr +m ) = ^ T ^ Kr A^^ir +m )] • (33) 
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We see that for a long time limit (here large n) and small amplitudes, the initial density contribution q(t) and the 
anomaly contribution A n (r) cancel each other out for N = 1. It explains why the lowest resonance (called sometimes 
semi — resonance j6j) is suppressed. No particles are produced inside a cavity. For higher resonance frequencies 
N > 1, the anomalous mechanism surpasses resonant enhancement of the negative energy of the initial state. We 
obtain again the exponential growth of the total energy and the formation of travelling narrow wave packets. 

The first exact analytical solution corresponding to a vibrating cavity system was presented in |lfl | for the resonance 
channel N = 2 and generalized in |16| for higher resonances. The solutions correspond to a family of mirror trajectories 
given by: 



L(t) = L + 



1 

UJ N 



. lo n AL 
sin — - — cos [u>Nt) 



(34) 



The picture of resonantly enhanced radiation is similar there except of the fact that the energy happens to grow 
quadratically with time. Such a power-like growth of the total energy is observed at the boundary of the frequency 
band Eq. I|23|) . For trajectories Eq. I)34|l . the billiard function can be derived exactly: 



/(t) = — arccot 

UJ N 



cot 



L) 



2 tan 



uj^AL 



L 



(35) 



where the branch of the multivalued function arccot should be properly chosen to have the billiard function increasing. 
It is straightforward to calculate the cumulative Doppler factor and the cumulative anomaly contribution: 



D n (r) = 1 + 2n 2 tan 



2 U N 



AL 



1 — cos (lu n (t + L))] + 2ntan UJn ^^ sm (cj n (t + L)) , 



(36) 



A„(r) = ^[l-D- 2 (r)] . (37) 

There are N periodic particle trajectories corresponding to starting points TQ m = (-N + 2m + 1)L/N where rn = 
0, 1, N — 1. But all D n (TQ m ) — 1. The periodic trajectories are neither positive nor negative. There is no Doppler 
factor since the reflections occur at the turning points where the mirror stops. We obtain the corresponding energy 
density: 

/ \ N 2 tt (N 2 — l)n { 9 „ WiyAL r .fj , Nl , >w ujmAL , . 1 ~ 2 , , 

B{ ) = ~ m? + 48L 2 I + ^ ~ C ° S {UJNT) ^ ~ 2n(_1) tan Sin ^ NT) ) (38) 

In this case, it is easy to show that the total energy grows quadratically with time, and the heights and widths of 
travelling wave packets are proportional to r 4 and r -2 respectively. 



IV. SUMMARY 



In summary, we have studied the resonance behavior of the electromagnetic field inside a vibrating one-dimensional 
(linear) cavity. Our approach exploits fully the analogy with the classical billiard of massless particles. The parametric 
resonance is determined by the existence of periodic particle trajectories. The billiard function can be used to establish 
recursive relations for all physical quantities. The resonant exponential enhancement and the concentration of energy 
into narrow wave packets is described in both classical and quantum field theories of cavities referring to the same 
aggregated Doppler factor. It is enough to calculate this factor, even numerically, in some finite interval of the length 
related to the period of cavity oscillations. A purely quantum phenomenon of particle production from the vacuum 
is described by the anomalous contribution in recursive relations. For higher resonances, this contribution dominates 
and the evolution of the system is insensitive to the initial state of the quantum field inside a cavity. However, for 
higher resonances the band structure defined by Eq. (|23|l is squeezed. From experimental point of view, the inequality 
Eq. H23fl challenges someone to secure either fine tuning or big enough amplitude with gigahertz frequencies. The 
resonance amplification and photon production from thermal field fluctuations p"?l ]. being practically more feasible, 
can be studied in the same way as well. 

The opto-mechanical resonance in vibrating cavities (dynamical Casimir effect) is the subject of numerous studies. 
The analysis of real, three-dimensional models shows that each photon mode can be truly described by some one- 
dimensional model. We have presented a new and simple method to analyze resonance solutions in such fundamental 
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models. Our new results: the mechanism of particle production and the mechanism of resonance amplification of ra- 
diation are separated and clarified, exact formulas for off-resonant solutions are given, the band structure is described, 
the puzzle with the lack of the amplification for the lowest resonant frequency is solved, it is also explained when the 
instability of the system (the energy growth) is either exponential or power-like. Moreover, the correspondence with 
classical optics and classical mechanics is presented here from a new point of view as far as the quantum parametric 
resonance is concerned. 
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